NASA/TM— 2001-210890 



Preliminary Study of Electron Emission 
for Use in the PIC Portion of MAFIA 

Jon C. Freeman 

Glenn Research Center, Cleveland, Ohio 


June 2001 


The NASA STI Program Office ... in Profile 


Since its founding, NASA has been dedicated to 
the advancement of aeronautics and space 
science. The NASA Scientific and Technical 
Information (STI) Program Office plays a key part 
in helping NASA maintain this important role. 

The NASA STI Program Office is operated by 
Langley Research Center, the Lead Center for 
NASA's scientific and technical information. The 
NASA STI Program Office provides access to the 
NASA STI Database, the largest collection of 
aeronautical and space science STI in the world. 
The Program Office is also NASA's institutional 
mechanism for disseminating the results of its 
research and development activities. These results 
are published by NASA in the NASA STI Report 
Series, which includes the following report types: 

• TECHNICAL PUBLICATION. Reports of 
completed research or a major significant 
phase of research that present the results of 
NASA programs and include extensive data 
or theoretical analysis. Includes compilations 
of significant scientific and technical data and 
information deemed to be of continuing 
reference value. NASA's counterpart of peer- 
reviewed formal professional papers but 
has less stringent limitations on manuscript 
length and extent of graphic presentations. 

• TECHNICAL MEMORANDUM. Scientific 
and technical findings that are preliminary or 
of specialized interest, e.g., quick release 
reports, working papers, and bibliographies 
that contain minimal annotation. Does not 
contain extensive analysis. 

• CONTRACTOR REPORT. Scientific and 
technical findings by NASA-sponsored 
contractors and grantees. 


• CONFERENCE PUBLICATION. Collected 
papers from scientific and technical 
conferences, symposia, seminars, or other 
meetings sponsored or cosponsored by 
NASA. 

• SPECIAL PUBLICATION. Scientific, 
technical, or historical information from 
NASA programs, projects, and missions, 
often concerned with subjects having 
substantial public interest, 

• TECHNICAL TRANSLATION. English- 
language translations of foreign scientific 
and technical material pertinent to NASA's 
mission. 

Specialized services that complement the STI 
Program Office's diverse offerings include 
creating custom thesauri, building customized 
data bases, organizing and publishing research 
results . . . even providing videos. 

For more information about the NASA STI 
Program Office, see the following: 

• Access the NASA STI Program Home Page 
at http:llivww.stunasa.gov 

• E-mail your question via the Internet to 
help@sti.nasa.gov 

• Fax your question to the NASA Access 
Help Desk at 301-621-0134 

• Telephone the NASA Access Help Desk at 
301-621-0390 

• Write to: 

NASA Access Help Desk 

NASA Center for AeroSpace Information 

7121 Standard Drive 

Hanover, MD 21076 


NASA/TM— 2001-210890 



Preliminary Study of Electron Emission 
for Use in the PIC Portion of MAFIA 

Jon C. Freeman 

Glenn Research Center, Cleveland, Ohio 


National Aeronautics and 
Space Administration 


Glenn Research Center 


June 2001 



This report contains preliminary 
findings, subject to revision as 
analysis proceeds. 


Available from 


NASA Center for Aerospace Information 
7121 Standard Drive 
Hanover, MD 21076 


National Technical Information Service 
5285 Port Royal Road 
Springfield, VA 22100 


Available electronically at 




Table of Contents 


Abstract 1 

Introduction 3 

Existing Electron Gun Software 6 

Analysis Summary 8 

Chapter 1 Emission Modeling 10 

Chapter 2 Low Voltage Diode 30 

Chapter 3 Particle Enumeration 41 

Beam Velocity Profile 52 

Conclusions 55 

Appendix 1 60 

References 63 

List of Symbols 65 


NASA/TM — 200 1-210890 


1 




Preliminary Study of Electron Emission for use in the PIC Portion of MAFIA 


Jon C. Freeman 

National Aeronautics and Space Administration 
Glenn Research Center 
Cleveland, Ohio 44135 


Abstract 

This memorandum summarizes a study undertaken to apply the program 
MAFIA to the modeling of an electron gun in a traveling wave tube (TWT). The 
basic problem is to emit particles from the cathode in the proper manner. The 
electrons are emitted with the classical Maxwell-Boltzmann (M-B) energy 
distribution; and for a small patch of emitting surface; the distribution with angle 
obeys Lambert’s law. This states that the current density drops off as the cosine of 
the angle from the normal. The motivation for the work is to extend the analysis 
beyond that which has been done using older codes. Some existing programs use 
the Child-Langmuir, or 3/2 power law, for the description of the gun. This means 
the current varies as the 3/2 power of the anode voltage. The proportionality 
constant is termed the perveance of the gun. This is limited, however, since the 
3/2 variation is only an approximation. Also, if the cathode is near saturation, the 
3/2 law definitely will not hold. In most of the older codes, the electron beam is 
decomposed into current tubes, which imply laminar flow in the beam; even 
though experiments show the flow to be turbulent. Also, the proper inclusion of 
noise in the beam is not possible. These older methods of calculation do, however, 
give reasonable values for parameters of the electron beam and the overall gun, 
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and these values will be used as the starting point for a more precise particle-in-cell 
(PIC) calculation. 

To minimize the time needed for a given computer run, all beams will use 
the same number of particles in a simulation. This is accomplished by varying the 
mass and charge of the emitted particles (macroparticles) in a certain manner, to be 
consistent with the desired beam current. 

Chapter 1 , Emission Modeling , gives a general overview of how an emission 
calculation should proceed, and develops equations that adhere to the basic 
constraints of the emission process. It ends by stating the number of particles for 
a simulation using this technique would be prohibitive. The remaining sections 
show alternate paths that will be pursued, to obtain outputs that will be similar to 
published experimental reports. Chapter 2, Low Voltage Diode, gives the 
analytical result for a 1 -dimensional planar diode with only 1.172 volts applied to 
the anode. This example demonstrates the magnitudes of the numbers involved. 
Then Chapter 3, Particle Enumeration , shows how to modify the particle mass 
and charge to reduce the run time to a reasonable value. Finally, the probability 
density for a reasonable model of an electron beam in other parts of the TWT is 
developed for analysis purposes. Appendix I gives a proof of constraints used in 
earlier sections. 
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Introduction 


This memorandum summarizes the study that was undertaken to determine 
if the electromagnetic solver, MAFIA, could be used to model the electron gun of 
a traveling wave tube. The preliminary conclusion is that it can be used; and the 
outline of how it will be executed is explained herein. 

A critical element of an electron gun is the thermionic cathode. The 
emission properties of any small portion of it may be viewed as the emitting 
surface of a 1 -dimensional diode. The solution for this 1 -dimensional diode 
model may be found in many references; and it is called the 3/2 power, or Child- 
Langmuir law. The simple solution is an elegant display of solving Poisson’s 
equation, but it does have some assumptions that make it difficult to implement 
in a program that uses a particle- in-cell (PIC) approach. 

To obtain a clearer understanding of the modeling difficulties, the actual 
process of the emission details should be considered first. The physical 
characteristics of the emission may be thought of in the following manner. At the 
specified cathode temperature, the electrons are emitted with a classic Maxwell- 
Boltzmann (M-B) energy distribution. As the anode potential is increased, the 
collected current increases and will eventually saturate. This condition is denoted 
as the “thermal limited regime”; See the figure below. 
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Figure 1 - (a) The J-V family of curves showing the space-charge-limited and the 
thermal-limited regimes. The parameter T is the Cathode temperature; Tj > T 0 . 
(b) Schematic of the 1 -dimensional diode. Near the cathode is the dense space- 
charge layer with its resulting decelerating Held. 


As the anode potential is reduced, the collected current will drop, and the 
diode is in its 3/2 law regime, which is called the space charge limited region. In 
this case the reduced electric field at the cathode does not drain off the thermally 
emitted particles fast enough, and a space-charge layer forms adjacent to the 
cathode. The space-charge layer acts to repel some of the thermally emitted 
electrons back into the cathode surface. This has the effect of retarding the slower 
particles, and only allowing the more energetic ones to climb the potential hill due 
to the space charge layer, and then proceed on to the anode. Experimentally [1], 
[2] it is known that the energy distribution of the ensemble that overcomes the 
retarding space-charge field is still M-B, but with a lower effective temperature. 

The space-charge layer causes a retarding electric field at the cathode 
surface, which decays rapidly to zero as one moves toward the anode. Beyond 
this point, the field is reversed in direction and begins to increase, and serves to 
accelerate the electrons toward the anode. The plane where the field is zero is also 
the point of the potential minimum (the virtual cathode). It is at this plane where 
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the classic boundary conditions are applied. Physically, the electric field is zero 
there, and the potential is just a few tenths of a volt below zero, but the carrier 
velocity is not zero; it is of the order of the thermal speed. The electrons leave the 
cathode with velocities around 10 5 m/sec (their thermal speeds). Since they are 
generally accelerated up to speeds of 10 6 or 10 7 m/sec in typical applications, their 
initial speeds at the real and virtual cathodes may be neglected. The potential 
minimum is very close to the actual cathode, and it is only a few tenths of a volt 
below the grounded cathode; so at this plane, the electric field, potential, and 
velocity are set to zero. These boundary conditions yield a solution that is very 
close to measured results. One problem is the fact that the charge density must 
be infinite (since the velocity is zero) at the virtual cathode, since the current must 
be constant in a 1 -dimensional problem. In a real diode, the charge density near 
the cathode is extremely large, but its value is not of particular importance; only 
the current- voltage relationship is relevant in device use. 

To model a cathode with a PIC code, one could emit particles at the cathode 
with a M-B energy distribution and a net current equal to the measured saturated 
thermal current density. Then for any applied anode potential, the space-charge 
layer would form appropriately and the collected current should match that 
measured to within numerical and experimental error. Unfortunately this straight- 
forward recipe is impossible to implement. A large fraction of the thermally 
emitted electrons are repelled back into the surface and do not contribute to the 
collected current, and are wasteful of computer resources. Also, the maintenance 
of the space-charge layer would be an overbearing tax on the computation time. 
The total number of particles to use gets to be unreasonable, and the run time 
would be prohibitive. Finding a way, or several ways, to circumvent this dilemma 
constitutes the discussion of this memorandum. 
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Existing Electron Gun Software 


While an exhaustive study of all available electron gun simulators has not 
been attempted, one of the most popular ones has been researched in some depth. 
The program is called EGUN, and was developed by Hermannsfeldt [3]. The 
program uses a “starting surface” which is several mesh cells away from the actual 
cathode. Current filaments are formed on the surface, and their values are varied 
until the net cathode current and anode voltage satisfy the 3/2 power law. That 
is the current is proportional to the 3/2 power of the anode voltage. The 
proportionality constant is the perveance of the gun, which is a function of the 
geometry of the electrodes. Notice this approach completely sidesteps the problem 
of forming the space-charge layer. In a real electron gun the electric field across 
the cathode surface varies due to the electrode geometry and applied potentials. 
This results in a variation of the space charge limited current across the surface, 
which is called variable cathode loading. This variation may change by 20%, with 
the largest emission density near the edge of the cathode, and the lowest emission 
density near the center. The program EGUN does calculate this variable loading 
phenomena. The electric field and potential are held to zero on the starting 
surface, and the iteration proceeds. The current filaments act on one another in 
a reasonable approximation, and magnetic focusing is incorporated. One level 
of approximation used on the initial EGUN calculation is the following. The 
documentation of the program dated in November 1979 (p. 85), says the program 
has difficulty determining the space charge on the beam axis. The calculated value 
is multiplied by the empirical factor of 5.5 to agree with some measured results 
(the details were not given). Notice the beam is “cold” (no thermal velocity), and 
thus noiseless. To simulate the noise in the beam, each filament is then subdivided 
into rays of varying strength and they are emitted at specified angles from the 
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normal to the starting surface. These filaments are quickly curved back toward 
the direction of the normal, as no crossing of the filaments is allowed in the 
calculation. It has also been stated [4] that the program determines the perveance 
to within a few percent, but has a somewhat harder time predicting the beam 
diameter. Unfortunately, for our analysis effort, the beam diameter determination 
is critical. This basic scheme is used by others [5], and is also incorporated in 
MAFIA as a macro entitled SCLE (Space Charge Limited Emission). 
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Analysis Summary 

With the physical concepts of the emission process in mind, the approach to 
model an electron gun using the program MAFIA is outlined here. The first step is 
to use the SCLE module to get the appropriate variation of cathode loading for the 
chosen electrode configuration. Next the cathode will be partitioned into segments 
(patches) where each segment is assumed to emit a uniform current density. Using 
this uniform current density in the patch, we will decompose it into the proper 
surface carrier density and velocity distribution to model a M-B distribution. 

Notice this step is emitting only those particles that have overcome the potential 
hill, (virtual cathode) and will be collected at the anode. The formation of the 
M-B distribution will be developed in the very long chapter entitled “Emission 
Modeling”. 

MAFIA requires the specification (in a patch) of the number of emitted 
particles, their charge, mass, energy (speed), position, and emission directions; 
this must be kept in mind in the “Emission Modeling” section. Notice the cathode 
is at temperature T } , with the appropriate M-B distribution for that temperature. 

At the voltage minimum (virtual cathode, and the location of the starting surface) 
the effective temperature is T eff- We must determine T eff, and the correct M-B 
distribution at this plane. A procedure is developed to provide this information. 

The emitted particles (in a patch) are grouped into 14 discrete energy bins, 
and the variation in the number emitted as a function of angle from the normal is 
determined exactly. After much trial and error it was decided to always emit the 
same number of particles for every simulation. The variation in current for 
different conditions is handled by varying the charge of the emitted particles 
(called macroparticles). In a M-B distribution the number of particles in a given 
energy (speed) bin is dependent on that energy. By altering the charge and mass 
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of the particles, we may use the same number of particles in each bin. The 
adjustments conserve both energy and momentum, and the change in mean-free- 
path is not severe. For example, we consider a low voltage diode later to get a 
feel for the magnitude of numbers involved in a calculation. If the diode is solved 
using electrons, their mean free path is about 5.89 microns. When we use our 
derived macroparticles, the path length is 15.67 microns. This variation will be 
assumed to be acceptable. The motivation for this plan is to keep each run to about 
the same time duration. From previous work it is known that using about 300,000 
particles per run requires about 4-6 hours of computing time. 

To verify the outlined scheme, a series of MAFIA computations for a known 
solution for a 1 -dimensional diode will be carried out. The details of the known 
solution are in the chapter “Low Voltage Diode”. 
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Chapter 1 Emission Modeling 

We know the energy distribution for the emitted particles is M-B both at 
the metallic cathode, and more importantly, at the virtual cathode (plane of the 
potential minimum). The program MAFIA requires the specification of: a) area 
of the emitting surface, b) the number of particles and their positions, c) the 
particle speed, and d) the particle direction. We may also specify the time step 
for successive bunch emissions. Given the temperature of a cathode 7), and the 
material, we can use Figure 2.3-1 of Gewartowski [6] to determine the saturation 
emission current density J 0 (Amps/cm ). The value J 0 is the limit in the thermal 
limited regime, (see [6] pages 42, 59, and 615). The sketch below indicates the 
typical J-V plot for a thermionic cathode. 




Figure 2 - The J-V curve for a single temperature, along with the Jo(6) for a small emitting 
patch. 

If we choose T } = T = 1100° K, we can assume the typical value of J Q = 0.5 A/cm . 
In the above sketch we have indicated an emitting patch and the current density 
in the direction 6 from the normal. Let n p ( 0, (p) and v (6, (p) be the number of 
particles emitted and their velocities in the ( 6, (p) direction, respectively. Here 6 
is the polar angle from the normal, and (p is the azimuthal measure. The average 
velocity, <v>, of an emitted particle (averaged over all directions) may be obtained 
from 
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( 1 - 1 ) 


^m(v} 2 =2 W t 

where m is the mass of a particle, and 2 W T is their thermal energy. W T is 

T 


W T =^r = 


■eV 


\q\ 1 1,600 

where k is Boltzmann's constant. For our example, 

2W t = 2kT = 2(l .38 x 1 0- 23 J /o K ){\ 1 00° K ) 
= 3.036 x\0- m Joule 
3.036x10 20 7 


. = 0. 1 895e V W = .095c V 

1.602x1 O' 19 J/, 


/eV 


Then 


W= 


j4W r _ 

" 6.072xl0" 2 °7 ’ 

V m 

9.108xl0~ 31 AT# 


-i>i 


Hence 


(v) = 2.582 xl0 5 m/ sec 

We arbitrarily choose l A (W T ) for energy of the slowest particle and 5.75 (W T ) for 
the energy of the fastest particle; then 

v sww = 1 -29 1 x 1 0 5 m / sec 
v fast = 6.191xl0 5 m/sec 

so we have a sense of the order of magnitudes with which we are dealing. Observe 
the ratio of fastest to slowest is only about 4.8. 

Assume no 0 dependence, then 

/„(0) = <?I«,(0)v(0) (1-2) 

where q = the particle charge. 

From reference [7] we have: 
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N r = total number of ions emitted per unit area per unit of time in the joint interval 
du dd d(j) is 

9 /j " fj 7 ~ , 

N r =^^-ide- hm ' r sinflcosfl du dd dtp (1-3) 

71 

where m = mass of ion (here just an electron), h = l/2kT (see pg 156), and w is a 
particular velocity. Note that the average energy of an emitted particle is 2kT 
(see eq. 1-1). This is in contrast to the value of 3kT/2 for a particle in the classical 
M-B gas. This is due to the fact that the group that has escaped, must have enough 
energy to overcome the retarding binding force at the cathode surface. For the 
particles remaining in the metal, their average energy is 3 kT/2. See Appendix I for 
the derivation of this result. Thus, the mean kinetic energy ( KE) in the metal may 
be expressed as 3/(4h). Integrate over 0; 

N r (u,6)= 4h 2 m 2 u 3 e~ hmu sinOcosOdudd (1-4) 

Observe this is eq. (2.4-8) of Gewartowski [6]. In eq. (1-4 ), u ranges from 0 to °°, 
but we know the major contribution is obtained during the interval 

V < u < V 
y SLOW “ w — * FAST 

Consider a group of particles with the same velocity (speed, call it u)\ then 
their distribution with 6 is 

N R (u l ,0)= (const.)- sinflcos# du'dd (1-5) 

Then the number per unit solid angle (U.S.A.) is 

NM . e )u.s.A . = const ■- — . du dO ( 1 - 6 ) 

Insmddd 

N l R (u 1 , 6 ) c , s A = const 1 • cos0 du 1 (1 -7) 

We interpret eq. (1-7) as stating "for a given particle speed, those with that 
speed are distributed as cos0." The units of N R (u\9) are [particles/unit 
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area]/[unit of time/U.S.A.] which equals [particle flux/U.S.A.]. Equation (1-7) is 
similar to the result of Gewartowski as follows. Start with his eq. (2.4-8) 


dP(u,6) = sin 0 cos 0 


^ m ^ 
kT 


u 3 e” U/ * T dud6 


d-8) 


where dP(u,d) is the probability that the emission velocity lies in the range u to u + 
du and makes an angle with the normal in the range 0 to 0 + dd. 

Integrate over all velocities (then the velocity dependence is removed). 

[dP(u,6)du = sin0cos0(c)~ j^u 3 e s du 

1 


= sin 0 cos 0(c)' 


= 2sin0cos 6d0 


. 2fe) = J 

dP(6 ) = 2 sin 0 cosd dd 

Which is eq. (2.4-1 1) of Gewartowski. Now we say 

J( 0) = [total emission current averaged over all directions] x 


d-9) 


dPj 0) 
U.S.A. 


dP{0) _ j 2sin0cos 0d0 
~ ° USA. ~ 0 27tsin0d0 

= — cost? emitted current density/U.S.A., which is Gewartowski's eq. (2.4- 
n 

12 ). 

For a given T, we know the range of emission energies. The fastest .5% of 
particles will have energy about 6W T . The distribution of particles as a function of 
energy is 

2N l~E~ s/ 

dN . = ——M = —e v "dE (1-10) 

E T 4n^E T 
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where E T = W T and E=W. Notice the symbol for energy will be either W or E\ we 
do this to conform to whichever author we are paralleling at a certain point in the 
analysis. Let p n = dN n / dp where p = E / E T . 

N TO t = particle density (particles/m 3 ). We will subdivide the range of energy into 
B energy bins. For now B is arbitrary. The velocity of a particle in binj is 


v ; = 


2 E 


m 


(Ml) 


Let fij be the particle density of those with energy E r The current from a patch is 


T = t A 

1 PATCH u o™ PATCH 


j 2 

where A PATC h is the patch area. Write J 0 in Amps/rtf and A PATC h in m 

I PATCH — J PATCH Q^PATCH (^ 1^1 ^ 2^2 ^ B ) 


d-12) 


(M3) 


where rtj = density of particles in binj j = 1 ,2. . .,B 

Vj = velocity of particles in binj 
We will show 


Yi zzoc AJ 

fl \ j i y T QJ 

Yi —oc AJ 

* 1 2 2 A T TOT 


(1-1 3a) 


YI =oc yV 

I^ B B 1 Y TQT 


Where the ©c are known constants. 

Then 

I PATCH ~~ PATCH ^ TOT 1 ^1 ^ 2 ^2 ^ ^ B ) 

but the Vj and are known, so 


( 


N. 


TOT 


particles 


m 


PATCH 


d^ PATCH ^ 00 j V j 


7 = 1,-, B 


(1-14) 


Once N TO t is found, rij are computed from eq. (1-1 3a). 
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From Millman and Seely [8], their section 5-15, we have 


KF = 2KF 

NORMAL TANGENTIAL 


(1-15) 


which states the KE normal to the surface is twice the tangential value, when 
calculated over the ensemble. See Appendix I for its verification. Let rip represent 
the particle density in energy birij and emitted at angle 0,. Let v y be the velocity of 
particles in buy. For calculation purposes we will choose 4 energy bins (j = 1 to 4) 
and 5 (/ = 0 to 4) emission angles. The sketch below illustrates the model. 


0 ° 


e, 


. 0 , 


/ 


/ 


/ / 


/ 

/ 


0 , 


0 . 


Figure 3 - The case for emission at 5 specific angles from a small patch. Emission occurs at 
0° as well as 0j through 0 4 . 


Eq. (1-15) becomes 

^{(«. o v i + « 2 o v 2 + « 3 o v? + n w v] )+ (n u v; + n 2] v; + n n v] + n 41 v 2 )cos 2 0, + 

(/ij, Vj 2 + n 12 y\ + « 32 v 2 + w 42 v 4 )cos 2 0 2 + (« 13 v 2 + ;t 23 v 2 + m 33 i> 3 + h 43 v 4 )cos : 0 3 + 
(n 14 v 2 + n 24 v; + « 34 v 3 2 + n M v] )cos 2 0 4 }= RHS 

Where RHS means the right hand side of eq. (1-15). 

But we know the velocities at the various angles are related by constants, i.e. 

V 2 =K 2 V v V 3 = K 3 V lt V 4 =RT 4 V, 
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n u =n l0 cos6 ] 

77 ,, ~ n io cosd 2 

n l3 = n w cos0 3 

n i4 = «,0 C0S ^4 
«JI = «30 C0S ^1 

« 32 = « 30 COS 0 2 

«33 = «30 COS ^3 
«34 = «30 COS0 4 


77,, — 77 20 COS0, 
77 , 2 = 77 , 0 COS0, 
77 2j = 77, 0 COS0, 

77,4 = «2o C0S #4 

77 41 = 77 ^ COS0, 

n 42 = U 40 COSd 2 
n 43 =n 4O cos0 3 

«44 = « 40 COS ^4 


So 


1 


{«10 V 1 2 + «20 (*T K + «30 (*7 V. 2 + «40 fc X 2 


+ (w, 0 v, 2 + ^20 K l v i + n 30 K 2 3 v; +n 40 K;v;)cos 3 9 l 
+ (n l 0 v x + n 20 K 2 v 2 + 77 30 7C 3 2 V 1 2 + n m K;v] )cos 3 d 2 + ( )cos 3 0 3 

+ ( )cos 3 0 4 }= v 2 X {l -I- cos 3 0, + cos 3 d 2 + cos 3 0 3 + cos 3 0 4 } 


( 1 - 16 ) 


Where 


X = 77, 0 + 77, 0 K 2 + 77 30 K 3 + U^k] (1-1 6a) 

Now the right hand side ( RHS) of eq. (1-15) is 

{( 77 , ,V 2 + 77,, V, 2 + 77 31 V 3 + 77 41 V 2 )sin 2 0, + (t 7 12 V 2 + 77 22 V 2 + 77 3 ,V 3 2 + 77 4 ,V 4 )sin 2 0 2 
+ ( 77,3 V , 2 + 77 23 v, 2 + 77 33 v 2 + n 43 v; )sin 2 6 } + ( 77 , 4 v 2 + 77 24 v 2 + 77 34 v 3 + n^v; )sin 2 0 4 } 

Or 

{(77 10 V 2 + 77, 0 7C 2 V, 2 + 77 30 7C 3 V 2 + 77 40 K‘ 4 V’ 2 )cOS0, sin 2 9, + 

( )cos0, sin 2 0 2 + ( )cos0 3 sin 2 0 3 + ( )cos0 4 sin 2 0 4 } 

= v, 2 { X (cos0, sin 2 0, + cos0 2 sin 2 0 2 + cos0 3 sin 2 0 3 + cos0 4 sin 2 0 4 )} 

Then 
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|v 2 X {l + cos 1 0, + cos 3 d 2 + cos 1 0 3 + cos 1 0 4 }= v 2 x(cos0, sin 2 0, 
+ cos0 2 sin 2 0, + cos0 3 sin 2 0, + cos0 4 sin 2 0 4 ) 

Then 

l + cos 3 0, + cos 3 0, + cos 3 0 3 + cos 3 0 4 = 2cos0,sin 2 0 1 + ■■ ■• 
Or 

• ^ + (fcos 3 0, -cos0, sin 2 0,)+ (4cos 3 0 2 -cos0 2 sin 2 0 2 ) 

+ (tCos 3 0 3 -cos 0 3 sin 2 0,)+(3cos 3 0 4 -cos0 4 sin 2 0 4 )=O 

Now define 

/(0) = 4cos 3 0 -cos0sin 2 0 

which is sketched below. 



d-17) 

d- 18 ) 
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Rearrange eq. (1-17) as 


/ (0, ) + / (0 a ) + / (0 3 ) + / 04 ) = "± C I ■ - 1 1 ' 7a) 


So we must choose the 9 l ,9 2 ,9 i ,9 4 to satisfy eq. ( 1 -1 7a) 


9° 

m 

9° 

m 

9° 

M 

9° 

f<8) 

1 

.499 

22 

.268 

43 

-.145 

64 

-.312 

2 

.498 

23 

.249 

44 

-.161 

65 

-.309 

3 

.495 

24 

.230 

45 

-.177 

66 

-.306 

4 

.492 

25 

.210 

46 

-.192 

67 

-.301 

5 

.487 

26 

.190 

47 

-.206 

68 

-.296 

6 

.481 

27 

.170 

48 

-.220 

69 

-.289 

7 

.474 

28 

.150 

49 

-.232 

70 

-.282 

8 

.466 

29 

.129 

50 

-.244 

71 

-.274 

9 

.458 

30 

.108 

51 

-.255 

72 

-.265 

10 

.448 

31 

.088 

52 

-.266 

73 

-.255 

11 

.437 

32 

.067 

j 

53 

-.275 

74 

-.244 

12 

.426 

33 

_ .046 

54 

-.283 

75 

-.233 

13 

.413 

34 

.026 

55 

-.291 

76 

-.221 

14 

.400 

35 

.0053 

56 

-.297 

77 

-.208 

15 

.386 

36 

-.015 

57 

-.302 

78 

-.194 

16 

.371 

37 

-.035 

58 

-.307 

79 

O 

00 

17 

.356 

38 

-.054 

59 

-.310 

80 

-.166 

18 

.339 

39 

-.073 

60 

-.312 

81 

-.151 

19 

.322 

40 

-.092 

61 

-.314 

82 

-.135 

20 

.305 

41 

-.110 

62 

-.314 

83 

-.119 

“ 21 

.287 ” 

42 

-.128 

63 

-.314 

84 

-.103 


NASA/TM — 2001-210890 


18 







85 

-.086 

87 

-.052 

89 

-.017 



86 

-.069 

88 

-.035 

90 

0.0 




The following 17 sets of angles that satisfy eq. (1-1 7a) will be listed for 
reference purposes. Notice that every angle from 1 to 90 is used. The sets were 
chosen randomly. 

Set Angles 

1 10,20,30,40,50,52,60,70,80 

2 5, 1 5,25,40,45,50,55,60,70,8 1 ,88 

3 15,20,30,43,48,53,58,68,87 

4 7,14,28,35,36,37,40,60,65,70,75,80,85 

5 2,12,31,42,47,51,56,62,74,86 

6 11,19,63,67,73,77,79 

7 3,8,54,57,63,67,72 

8 1,59,66,71,84 

9 4,61,73,76,77 

10 23,64,69,82,89 

1 1 22,24,26,41,44,46,49,78,79,83 

12 6,18,34,38,39,45,46,57,61,75 

13 9,52,64,73,77,85,88 

14 13,32,68,71,73,81 

15 16,33,48,64,74,84,88 

16 17,29,42,62,68,76,89 

17 21,27,34,44,56,67,77 

Some 4 angle sets are: (notice their sum should be -0.5) 

15,60,66,71 their sum = -.506 
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20,48,57,70 their sum = -.499 
25,62,75,80 their sum = -.503 
30,68,77,84 their sum = -.499 
Recall 

(«,o +«u +n l2 +n 13 +n l4 )=n l 
Or 

« 10 (l + COS0, + cos0 2 + cos0 3 + cos(? 4 )— /i, 

n, (1-19) 

1 + , 

And n 2 o , n 30 , n 40 can be expressed similarly. Thus all 16 bunches 

riji, i,j— 1 ,2,3,4 are determined. For each angle set, we will distribute them 

uniformly in phi (the azimuthal angle). 

Now we calculate the rij. Recall from eq. (1-1 3a) 

fly 00 j N j 0T 

flj -'x-j N tot 

The energy distribution function is shown in the plot below 



Figure 5 - The Maxwell-Boltzmann energy distribution function plotted versus the 
normalized variable r)=E/E T . 
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Consider a small patch which emits N T ot particles over some specified 
time interval. Subdivide the emitted particles into 14 groups (now we will use 
14 energy bins), with the particles in each bin having the same energy, (and thus 
speed). By graphically making the partitions, we may calculate the number of 
particles in birij 

0 - 20 ) 

Evaluating eq. (1-20) in general terms gives 

n = T ^ e "’ dr l 

Where we have suppressed the coefficient 2N TOT / 4tz 
Let 


So 


Let 


Or 


Then 



* = rj. 


2 xdx = drj 


n — 2 1' 2 x 1 e dx 


( 1 - 21 ) 


f = e~" 

f"=-2xf- 2/ = 4x'-f -2/ 

*7 = i(/"+2/) 
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2 \ x * x 2 e X dx 


= i/‘ 


+ 


J 


=i|;(/"+2/)fa 

= ■([/" (.t, ) - / 1 (a;, )]+ J:" <T’’ * 



e ' dx = i[2xj(x,)- 2x,f(x 2 )]-i 


4n 

~Y 


[erf(x 2 )-eif( *,)] 


So 

fZ 

n - xj (a-,) - x 2 /(x 2 ) + [e>/(x 2 ) - ^//(x,)] 

Thus 

n, =l^n = N nr lj={x l e-’ ; -x 2 e- l )+erf(x 1 )-erf(x l ) J 


(1-22) 
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Our values are: 


1 1 

X 


xe ' 

erf(x) 

0 

0 

1 

0 

0 

.5 

.70711 

.60653 

.42888 

.68260 

1 

1 

.36788 

.36788 

.84270 

1.5 

1.2247 

.22313 

.27327 

.91679 

2.0 

1.4142 

.13534 

.1914 

.95449 

2.5 

1.5811 

.082085 

.12978 

.97464 

3.0 

1.7321 

.049787 

.086236 

.98569 

3.5 

1.8708 

.030197 

.056493 

.99185 

4.0 

2.0 

.018316 

.036632 

.99532 

4.5 

2.1213 

.011109 

.023566 

.9973 

5.0 

2.2361 

.0067379 

.015067 

.99843 

5.5 

2.3452 

.0040868 

.0095844 

.99909 

6.0 

2.4495 

.0024788 

.0060718 

.99947 

6.5 

2.5495 

.0015034 

.0038329 

.99969 

7.0 

2.6458 

.00091188 

.0024127 

.99982 


Then 
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+ .6826 j = . 1 9866 N Tor 

'h = N rm { }= .2289 N tot 

n i = . 1 8085iV rO7 . 

n 4 = . 1 3008A^ ror 

7z, = .089682iV ro7 . 

n 6 = .060185A ror 

n, = .039722 N T0T 

n s =.02588 IN tot 

n q = .016724A ror 

n l0 = .0\012 N tot 

n n =6.8466x1 O' 1 N T0T 

7i 12 = 4.3436 xlO'Xor 

72,3 =2.7464 xlO" 3 A ro r 
72 14 = 1.7326 x\0~ y N TOT 


«,=A' ror {^("- 4 2 88 ) 


As a check, 

S ", =-9971N ror 

Thus the final table is: 


(1-23) 


Bin 

*7, 

Oj/NjOT 

1 

.25 

.19866 

2 

.75 

.22893 

3 

1.25 

.18085 

4 

1.75 

.13008 

5 

2.25 

.089682 
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6 

2.75 

.060185 

7 

3.25 

.039722 

8 

3.75 

.025881 

9 

4.25 

.016724 

10 

4.75 

.01072 

11 

5.25 

6.8466x10 3 

12 

5.75 

4.343x1 0’ 3 

13 

6.25 

2.7464x1 O' 3 

14 

6.75 

1.7326x1 O' 3 


Where rj, is the normalized energy parameter at the center of the partition. 
The emission velocity is developed as follows. 

Imv 2 = qE . ; write E j =2E T 7j i 


Where E T = ^qq ; the factor of 2 is to account for those that have escaped. 


■kmv 2 . = 2 Tj. 


f Y \ 
11,600 


= 27 ijkT 


(1-24) 


or 


\4kT- 
I — n j = 

m 


-23x1^ 


4(1.38x10 ) 
9.108xl0 -31 

= 7.785x10' yj^jT ml sec 




(1-25) 


For our case, T =1 100° K, thus v. = 2.582x10 ^jrj. 


d-26) 
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Bin 


Oj/Ntot 

Vj (m/sec) 

1 

.25 

.19866 

1.29xl0 5 

2 

.75 

.22893 

2.2361xl0 5 

3 

1.25 

.18085 

2.8868x1 0 5 

4 

1.75 

.13008 

3.4157xl0 5 

5 

2.25 

.089682 

3.873 “ 

6 

2.75 

.060185 

4.282 “ 

7 

3.25 

.039722 

4.655 “ 

8 

3.75 

.025881 

5.0x1 0 5 

9 

4.25 

.016724 

5.32 “ 

10 

4.75 

.01072 

5.627 “ 

11 

5.25 

6.8466x10 3 

5.92 “ 

12 

5.75 

4.3436xl0' 3 

6.19 “ 

13 

6.25 

2.7464x1 0' 3 

6.46 “ 

14 

6.75 

1.7326xl0' 3 

6.71xl0 5 


Note Vpast / v slow — 5.2 
Gewartowski shows 


7(e)=i cose A *” /cwl 

v ' n U.S.A 

Assume all particles are identical, then 

J(6)=q^L[nX e ’ < f>y i +« 2 (0,0V 2 + -+h 14 (0,0>' 14 ] 

Where W/ (0, 0) = the number of particles in direction (6, 0) with speed V/ 
Then 

11^(0, <p)d6d(f) = n l = k { N tot 


(1-27) 

(1-28) 


NASA/TM— 2001 -210890 


26 



In general 


ll n,(,e,<l>)ded<t> = n l =K l N T0T (1-29) 

It is shown in Richardson [7] that for a fixed velocity, the angular distribution is 
uniform in 0 and cosine in theta (9). Thus 

n j (9,(j)) = AjCOsO (1-30) 


Then 


k N = P : (2tt )4 cos 9 d9 = 2 kA 

J TOT J9= 0 v r J J 

K ; N t 


. a — 

1 2k 


Observe eq. (1-30) satisfies a fundamental constraint; namely 


KE =2 xKE, 

normal tan gential 


For a given velocity 


m t*A i \~ mA v T , 

KE norma ,=—[ A j cos9{v j cos9)d9=—^—f 0 cos 9d9 


=^(K) 


d-31) 


d-32) 


d-33) 


KE, 


mA.v 2 


_ J 


tan gemial 


L j^ 2 cos 9 s'm 2 9d9 = — ^ (1-34) 


Observe eq. (1- 33) is twice eq. (1- 34), which validates the procedures used. Now 
combining eqns. (1-27) and (1-28) 

J(9) = —cos9 = ^{(A. cos#)^ + (A 2 cos#)i\ H h (A l4 cos6>)v 14 } 

k 


Or 
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(1-35) 


{KjA^Vj + k 2 N IOT v 2 + • • 

n 2 n 

2 ;= i 

Now we can determine Njot given J G and T. 

A typical example might use J 0 = 5 Amp/cm 2 


14 

Evaluating ]Tk\v ; 


1 


K.v, = 2.565 xlO 4 

k 9 v 9 = 

K" 2 v 2 =5.1 19x10 4 

VlO 

k 3 v 3 = 5.221xl0 4 

K n v n z 

jc 4 v 4 = 4.443 xlO 4 

K 12 V 12 

K’jVj = 3.473 xlO 4 

^*13^13 ’ 

jc 6 v 6 = 2.577 xlO 4 

^14 V 14 

K 1 v 1 = 1.849 xlO 4 

14 

K’gVg =1.294xl0 4 



From eq. (1-35) 


N 


TOT 


2 J. 1 
q 


= 2.15x10" «/ 
/ m 
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' 4 " ^14 N T OT ^14 } 


= 5 x l(f Amp/m 2 

8.902x10* 

= 6.032x10* 

= 4.051x10* 

= 2.689x10* 

= 1.773x10* 

= 1.162x10* 

= 2.9 x 10 s 



4 = 6.798xl0 16 
4 = 7.834xl0 16 
4 =6.188xl0 16 
A 4 = 4.451 x 10 16 
4 = 3.069 xlO 16 
4 = 2.059 xlO 16 
4 =1.359xl0 16 


4 = 8.856xl0 15 
4 = 5.723 xlO 15 
4 0 =3.668xlO ls 
A„ =2.343 xlO 15 
4, =1.486 xlO 15 
4, =9.398xl0 14 
4 4 = 5.929 xlO 14 


While the above is without approximation (except for the finite number of 
partitions), the numbers are too large for MAFIA. There are many constraints on 
the emitted particles when using the program; some of them are: There is no easy 
way to incorporate the cosine distribution with 6 . Instead, we’ll use a parabolic 
distribution that closely matches a cosine. The number of particles/bin cannot be 
too large, or the run time will become prohibitive. The number of patches may be 
large, especially when modeling near the edge of the cathode. The question as to 
whether one should use systematic emission versus random emission is a very 
perplexing problem. The systematic case develops lower computation noise, but 
a random process seems to be what one needs. However, how does one produce 
such a condition, yet guarantee the emission is theoretically correct? Another 
constraint is the use of 4 to 5 particles per mesh cell to obtain good results. These 
along with other conditions will be explored and overcome in the succeeding work 
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Chapter 2 Low Voltage Diode 

To obtain a clearer perspective of the magnitudes of the parameters for a 
thermionic cathode, an analytical study of a typical 1 -dimensional diode was 
undertaken. The following example is close to that given by Beck [9], pg. 173 
(the example has some errors). The diode has plates with surface area of 0.2 cm" 
and separated the distance d, as shown below. 



Figure 6 - The geometry of the 1 -dimensional diode. 

The current drawn is 1 .0 mA, the current density is 5 x 10' 3 A/cm , and the thermal 
limited value J a is 5.0 A/cm 2 . At the cathode the parameter rf is 

Tl = ri c = In (J 0 / J) = In 1 0 3 = 6.9078 

From the table on page 527 of [9], we see this corresponds to the parameter 
-£ c = 2.509. Where £ is defined by 

£ =9.174x10 5 T^V7(x-xJ (2-1) 
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where x, x m are in cm, J is in A/cm 2 , and T is in Kelvins. The voltage minimum is 
located at x m . This equation is taken from the original paper by Langmuir [10]. At 
the cathode, £ = £ c . and x = 0. Then 

- 2.509 = -3.3963xl0 2 x„ 

m 


or 

x =7.4x1 0” 5 m = 0.074mm 

m 

Then the voltage at the minimum is 

V = In 1 0 3 = -0.6547 1 volts 

1 1,600 


At the anode 

£ a = 3 .396 x 1 0 2 (.05 - .0074) = 1 4.468 
From the table in [9] we read T] a = 19.27 
Now 

1 9.27 = ii^(F a +.65471) 

1100 

V = 1.172 volts 

a 

Where V a is the required anode voltage. 

With this information we may generate a plot of the potential through the 
diode versus x. The table is: 


x(cm) 

( 

n 

V (volts) 

0 

-2.509 

6.9078 

0 

.000625 

-2.301 

3.5 

-.323 

.00125 

-2.089 

2.286 

-.438 

.0025 

-1.6642 

1.12 

-.549 
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.005 

-.8151 

.2 

-.636 

.006 

-.4755 

.065 

-.649 

.007 

-.1359 

.0045 

-.6543 

.008 

+.204 

.001 

-.6546 

.009 

+.543 

.067 

-.648 

.010 

+.883 

.17 

-.6386 

.020 

+ 4.279 

2.8 

-.389 

.030 

+ 7.676 

7.2 

+.028 

.040 

+ 11.07 

12.7 

+.549 

.050 

+ 14.47 

19.27 

+ 1.172 


The plot of the voltage is shown below. 



Figure 7 - A sketch of the voltage V(jc) through the diode. The value at the minimum 
position (ar OT ) is about -0.655 volts. 
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We could numerically differentiate V(x) to obtain the electric field, E but to 
check for consistency between various treatments of the 1 -dimensional diode, [11], 
we will use an analytic formula for E from [9]. Forv <x m the field is given by 

E 1 = *:[<•"- 1 + e'erfjj - ifyif] 

For v > x m , we have 

T = K^' -l-e'erfjij + 2(y.f] 

Where E is in V/cm, K = 6250J -Jf = 1.0364 x 10 3 C/cm 2 . The table is: 


x{cm) 

n 

Jl 7 

er f(yfn) 

E(V/c m) 

0 

6.9078 

2.6283 

.9998 

1.4382xl0 3 

.000625 

3.5 

1.8708 

.9918 

255 

.00125 

2.286 

1.512 

.967 

131.3 

.0025 

1.12 

1.058 

.865 

60.4 

.005 

.2 

.447 

.466 

17.17 

.006 

.065 

.255 

.281 

9.07 

.007 

.0045 

.067 

.075 

10^ 

.010 

.17 

.4123 

.438 

-11.67 

.030 

7.2 

2.683 

.99985 

-48.06 

.040 

12.7 

3.56 

1 

-55.96 

.050 

19.27 

4.3898 

1 

-64 


The first and last E -field entries are those at the cathode E c and anode E a . The 
magnitude of the field at the anode is 64 V/cm , and from the graph of the potential, 
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the slope at the anode is about 1 volt over .015 cm, or 66.7 V/cm. The agreement 
is well within tolerance. 

Next we may determine the charge density p(x ) as follows. For 
convenience we first determine it in the accelerating part of the field; 

For x>x m 

r = 4?" - 1 - e'erfjrj + 2(?/)"] 

Let f l (r)) = e , erf{jrj) 


A'(v)= 


4n 


+ etfyfri e n =-= + e n erf^r\ 
yjn 


/(*?)= ^ -l-fXv)+ -J=ri* 

V7T 


f x {n)=e" -f;(n)+ 


2 1 1 


4n 2 Jrj 


- e 


i — 4 n + i — 

yin ^Jn K] 


e n ei 


But 


2E^=Kf(n)f 

ax ax 


dp dr\ dV — dr\ 
~ — —L - 


dx dV dx 


2E^=Kf‘(n 

dx 


-E 


dV 
drj 
dV 


2—--K 


4n 


-e +n erf4n + ~r 




drj 


dV 


where £ g is the permittivity of free space. 
Or 
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2 


p(x) = - 


2 

Ke 


44 

2 


-e^erf-Jr) + 


4™) 


dr\ 


44 


e e 


rf4n + 


444 


dV 

( 1 1 606 


1100 


Note 


Where 



r l 1606 ^ 
, 1100 , 


K = 1 .0364 x 1 O’ (y/ |=1.0364xl0 lV V 2 
V cm } / m 

p(x) = -4.8409 x 1 O'" |/' (rj)| cou f n 3 
% = 3.3962 xl0 2 (x-x m ) 


Our table is below. 


x( cm) 


V 

44 

erf 44 

|/'(77) 

p{dm') 

.0075 

.03396 

.00027 

1.64E-2 

— 

— 

~ 

.0076 

.0679 

.001 

3.16E-2 

.338 

17.38 

-84.1x10"* 

.008 

.2038 

.0096 

1 

9.8E-2 

— 

— 

-4.59x1 0 4 

.009 

.5434 


.259 

.284 

1.0993 


.010 

.883 

.17 

.412 

.438 

.906 


.020 

4.28 

2.8 

1.67 


.49522 

-2.4x10"* 

.030 

7.68 

7.2 

2.68 

.99985 

— 

— 

.040 

11.07 

12.7 

3.56 

1 

— 

— 

.045 

12.77 

15.7 

3.96 

1 

— 

— 

.050 

14.47 

19.27 

4.39 

1 

.99988 

-4.84x10"* 


N AS A/TM— 200 1-210890 


35 













Now for x < x m . 


So 


tsal 

l-J 

II 

e n - 1 + e n erf yjrj - 2 

'it 



, n ) _ 


p(x) = 




x(cm ) £ r\ fl 07 )| p(coul/m 3 ) 

0 -2.509 6.9078 2.001 x 10' .9687 


Now we estimate p from the slope of E . 


The final table is 



\ J 


x(cm ) 


3.125x10"* 

-.1673 

9.375x10"* 

-.01753 

1.875x10 3 

-.005 

3.75xl0‘ 3 

-1.531x10 3 

5.5x1 0' 3 

-7.1 7x10"* 

8x1 O' 3 

-4.59x10"* 

2x1 O' 5 

-1.61xl0 5 
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3.5xl0' 2 

-6.9947x1 0 5 

4.5x10“ 

-7.1 19x1 0' 5 


Now we determine the total space charge between the plates. 

coul/ _ 


d E / \ 

£ o—=p W 

dx 


m 


So 


e^-iy^dx 

- dx zi w , 

V 

coul/ 

/ tn 


e\E(x b )-E(x„)]= Q / 


Where x a and Xb are arbitrary points, and A = cross-sectional area. 
For x < x m 


E = K 


e n -l + e 




x = 0, E(0) = 1. 4382x1 0 s %,; jt = at„ E(xJ = 0 

So the charge between x = 0 and x = x m is 


-1.4382x 10 3 V xlO 2 — 1 = 

/cm m J 

A - .2cm 2 =2x10 ~ 5 m 2 

So £, = (2xl0' 5 m 2 )(8.854xl0~ 12 %)(- 1.4382X 10 5 %) 

= -2.547 xlO -11 coul => 1.59xl0 8 e~ 

For x > x w 


0, 


— 2 

£ = A' 


-1- e^erf /rj + 2 


'l±#' 

K 
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E - K\\ - 1 - 0 + 0] = 0 


At x = x m 


Hence E(x m ) = 0 as before. 

Aix = x a (the anode) E(x a ) = -64 



:.Q 2 = -.1133x10" coul => 7.08x10V" 
Then the total charge is 


Qtot =Q [ +<2: =-2.66x10" com/ => 1.66x10V 

We may now ascertain the velocity field: 

J = 50 Am P/ mi = p(xHx) 


'(*)= 


50 


pM 


m/sec 


The table is: 


x(cm) 

v(m/sec) 

3. 125x10"* 

2.99x1 0 2 

9.375xl0' 4 

2.85xl0 3 ” 

1.875x10"* 

9.956x1 0 3 

3.75x10 3 

3.27xl0 4 

5.5x10 3 

6.97xl0 4 

8x10 3 

1.089xl0 5 

2x1 0' 3 

3.1xl0 5 

3. 5x1 O' 2 

7. 148x1 0 5 

4.5xl0' 2 

7.02xl0 5 
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Note the thermal emission velocity is about 4.0E+5 ml sec. So the values in the 
above table are very reasonable. Notice the velocity varies by about a factor of 
1000. 

The magnitude of the electric field E(x) , across the diode is shown below. The 

cathode value is 1,438.2 volts/cm, while that at the anode is -64 volts/cm. The 
field passes through zero at 0.074 mm from the cathode. 



Figure 8 - A plot of \E\ through the diode. 
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Below is the charge density p(x ) and the velocity field v(x) across the diode. 



Figure 9 - The charge density p(x) and velocity field v(x) for the problem at hand. 


This completes the calculation of a 1 -dimensional diode in a regime where 
the 3/2 power law is not applicable. That regime is when J/J Q is between very 
roughly 0.01 to 0.5, rather than our case of 0.001. It is a useful example as the 
voltage minimum is away from the plane of the cathode and the charge density is 
not extremely large near the cathode. Here the density drops 3 orders of magnitude 
from the cathode to the voltage minimum (the virtual cathode), and the applied 
voltage is only 1.172 volts. We are drawing 1 mA, which is not too far from some 
actual TWTs (at least order of magnitude). Reference [1 1], by Kirstein, Kino, and 
Waters, has a good treatment of the problem, albeit in a notation that is unique. 

See refs. [12] through [15] for more information on this problem. 
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Chapter 3 Particle Enumeration 

We now attempt to determine the way to model the 1 -dimensional diode 
using a finite number of particles emitted over some specified time interval. First 
we divide the space between the plates into equal segments A d, and over the 
segments we use At = Ad / Av, where Av is the change in the velocity field over 
a segment. By using 9 such equal length segments we ascertain the transit time 
is about 27 nsec (for the Low Voltage Diode studied earlier). Gewartowski [6], 
pg. 629 gives the time of flight for the 3/2 law regime as 


<f = 


L nJ _ 


-i* 


which for our case is 1 .45 nsec. The difference of more than an order of 
magnitude shows how far we are from the normal regime. 

Recall, we would like to emit particles at the virtual cathode (plane of the 
voltage minimum). Help in this task is provided by [1 1] on pg.275, eq. (2.44). 
This expression gives the variation of potential with distance near the voltage 
minimum 


J in Amp/m 2 


V-V = 


4n\. 82x10" 


2e 4f 




dV 

dx 


4n_ 1.82XHT* 

2 £ 0 Jf 



(3-1) 


(3-2) 


The next derivative gives the charge density, which we see is not a function of 
position about the virtual cathode. 

d 2 V 4n 1.82x10^ T 

— - = = J i 

dx e dT 
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Observe eq. (3-2) shows that the electric field varies linearly through the virtual 
cathode plane. Then 


p(0= 


-^1.82x10 r4 J 


IJL 

1100 


-1.82x10 




-4.86xl0" 4 9 


m 


Notice the analytical result in the section “Low Voltage Diode” gave the charge 
density of -7 x 10 -4 coul/m J , so our results show remarkable consistency. Since 


p( x m )= -82 x 10- 4 [(- p(x m ))(v(x m ))] 
(v(x m )) = J^-5.49 x 1 O' ml sec 


(3-4) 


Which states the average velocity at x m is independent of all other parameters 
except the cathode temperature. For typical ranges of T we find 


T(K) 

v(x m ) m/sec 

1000 

9.79xl0 4 

1050 

l.OxlO 5 

1100 

1.03xl0 5 

1150 

1.05x10 s 

1200 

1.07xl0 5 


So <v (x m ) > is nearly 1 .0 x 10 5 m/sec for all temperatures and applied voltages, 
and currents. For reference, the velocity values in bins 1 and 14 are 1.29 x 10 5 and 
6.71 x 10 5 m/sec. 

Observe that the energy needed to overcome the retarding field of the space 
charge layer is 
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KE = q(0.655volts) = 1 .048 x 10' 9 Joules 
The velocity required is 


—mv 2 =1.048x10 19 
2 

.\v = 4.8xl0 5 m/sec 

So only those particles in bins 8- 1 4 can escape the region and proceed on to the 
collector. This amounts to only 0.07 of the total number emitted per unit time. 

Now consider a diode carrying the current i . Then the current in a given 
patch is 


= i 

p N At 

P 


(3-5) 


Where N p = the number of patches chosen for the entire cathode. Then 

Ae,=i,Af = Xn,Ag (3-6) 

Where n p = the number of particles per patch, and AQ is their net charge. Here 
a<2„ is the total charge emitted in a patch in the specified time interval At . 


Now 


A<2_ = A Q ] + A Q 2 + • • ■ + A Q l4 


where 


A Qj AQ I4 n l4 q l4 

Since there are 14 energy bins. Here qj is the charge value in a given bin. 
Consider biny, its total kinetic energy (KE) is 


(3-7) 


(3-7a) 


f 


KEj =», 


M 


i \ j 


(3-8) 
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Where nj, M jy v ; are the number of particles, their mass, and their common speed 
respectively. If the particles in all of the bins had the same mass m (this is the 
actual case in a real device), then 


KE l4 

KE j 



(3-9) 

(3-10) 


re-arrange 


KE. n. ( v, 3 


KEl 4 «.4 


v v - 




Now let us use n i4 particles in binj, but keep its KE the same. 


KE J = n l4 


{ -M 
2 J ' J 


We retain the original KE by assigning the particles a new mass Mj 
From eqns. (3-11) and (3-12) 


fJjKE l4 - K . - n l4 


- M .\ 2 
2 3 ) ! 


Or 


n 

V / 

j 

J 

«14 

v l{ 


m u V i4 n i4 


( 


= n 


14 


-M \ 2 

2 1 j 3 


(3-11) 


( 3 - 12 ) 


(3- 13a) 


(3- 13b) 


Or 


M =— m 14 (3-14) 

«,4 

Now the ratios n/n 14 are known from the M-B distribution, and for our chosen 
energy bins, they are: 
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(3-15) 


71, ln l4 = 1 14.66 
n 2 /n u =132.13 
n 3 / 72 14 =104.38 

72 4 / 72,4 = 75 

72 5 / 72, 4 =5 1 .76 

72 , / 72 .. = 34.74 


72 8 / 72, 4 = 14.94 
72y / 72,4 =9.65 

72.0 / 72 , 4 = 6.19 
72,,/ 72, 4 =3.95 

72.1 / 72, 4 =2.51 

72 n / 72 M =1-59 


72, / 72,4 = 22.93 

Since MAFIA keeps the charge to mass ratio constant (equal to that of an electron), 

< 7 , _ < 7,4 


M 


772 


14 


Hence 


M 


4, 


Returning to eq. ( 3-7) we observe 


722, 


72, 


r i 4 


(3-16) 


Afi,- — 


722, 




(3- 16a) 


And so forth; Thus 


AS = 


M, M, 


M 


\ 


13 


v 772,4 722,4 


722 


14 


+ 1 P, 472, 4 


(3-17) 


Where we have set the number of particles in each bin to be n 14 , as stated after eq. 
(3-1 1). Now from eq. (3-14) M / mu = n, / nu , so the term in parentheses in eq. (3-17) is 
known 


72, 


72 


A 


13 


+ 1 


^ 72,4 


72 


14 


= 575.43 


Where we have used eq. (3-15). Thus eq. (3-17) becomes 

A(2„ = 575.4327, 4 t2 , 4 =i At 


(3-18) 


(3-19) 
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Now eq. (3-19) is our basic working equation. The current from the patch i p 
is known, the time interval At may be chosen compatible with MAFIA’ s 
constraints on run time, etc. The number of carriers may be chosen, then the 
charge on the particles (macroparticles in MAFIA) is determined. Eq. (3-19) 
determines the charge on the fast particles, which are in bin 14. The mass and 
charge for the remaining 13 bins are found using eqns. (3-14), (3-15), and (3- 16a). 

Before q ]4 is determined we must choose n I4 . The total number of charges 
emitted during the interval At is 

N^=n„(l4)N p (3-20) 

Where N T otc is the total number emitted for all patches. For a given energy bin 
we know the particles are distributed as cos 6, so we choose the following recipe. 
Let 1 particle emerge at 85° and let the number increase as 6 — » 0 . The table is: 


e° 

# particles 

cos 6 

85 

1 

0.0872 

80 

2 

0.1737 

75 

3 

0.2588 

70 

4 

0.342 

65 

5 

0.4226 

60 

6 

0.500 

55 

7 

0.574 

50 

7 

0.643 

45 

8 

0.707 

40 

9 

0.766 

35 

9 

0.8192 
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30 

10 

0.866 

25 

10 

0.906 

20 

11 

0.9397 

15 

11 

0.9659 

10 

11 

0.9848 

5 

11 

0.99619 

0 

11 

1.000 


Then n ]4 = 136 particles/energy bin. The number of particles processed in a 
complete MAFIA run is 

W rOT _„„=(136)(14)W,^ (3-21) 

Where N cyc is the total number of emission cycle intervals. 

We can apply this to our low voltage diode problem as follows. The total 
charge between the cathode and the voltage minimum is 

Q ] =-2.5468x10 " C => 1.5917 xlOV 
While that between the virtual cathode and the anode is 
Q 2 =-.1133x10 " C=>7.0833xl0V 
Q tot =Q 1 +Q 2 =-2.66x 10-" C=>1.6625xlOV 

Since the time of flight is about 27 nsec, choose At = 2.7 nseconds. Choose 
10 patches, then i p = 10 -4 Amp, and from (3-19) 

AQ p = 10- 4 (2.7 x 10 -9 ) = 2.7 x lO'^C 

So the amount of charge developed/emission cycle is 

10A<2 p =2.7x1 0 12 C 

The total charge between the plates is then generated in 
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2.66X10" 11 


2.7 xlO” 12 

Also from (3-19), the charge q I4 is 


9.85 cycles 


< 7.4 


2.7x10‘ ,3 C 

(575.43X136) 


= 3.45xlO” I8 C => 21.56c” 


If we use 50 cycles, the total number of particles processed is 

N tot -run = (1 36K1 4)(1 0)(50) = 9.52 x 1 0 5 particles/ run 
So 9.52 x 10 5 particles emitted during 135 nsec. This value corresponds with that 
used by C. Kory and K. Vaden using the program. C. Kory uses about 2.8 x 10 5 
particles per run. 

At this point we have all the information to emit particles from the virtual 
cathode. We now must determine the correct M-B energy distribution at this 
plane. We know the saturation current is 10 3 times larger than that collected. Only 
10 3 of that emitted by the actual cathode is actually collected; the remainder is 
pushed back into the cathode, only to be re-emitted again. In other words, as we 
move from the cathode toward the voltage minimum, the carrier density decreases 
rapidly, but the imbalance of carriers moving in opposite directions still maintains 
the net 1 mA current toward the anode. At the virtual cathode and beyond, 
particles are only moving toward the anode and they constitute the 1 mA current 
at every plane. 

The velocity distribution for each bin was given in the Emission Modeling 
section, eq. (1-26) and subsequent table. Eq. (3-4) above shows the average 
velocity at the minimum to be 

{v(x m )) = J^5 .49 x 1 0 3 m / sec (3-22) 


And for our temperature, T = 1 100° K, we have 
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(v(jc m )) = 1 .03 x 1 0 5 m / sec 

From [8], the most probable speed occurs when r\ . = 1, (see eq. 8.32, pg 224). 
Then from eq. (1-26) of the Emission Modeling section we find the most probable 
speed = v r = 2.582 x 10 5 m/sec. Then the average speed v is 1.128 v c = 2.9125 x 
10 5 m/sec. The rms value is 1.224 v c = 3.16 x 10 5 m/sec. Also from eq. (1-26) 

v c = 7.785x1 O' Vr ( since77;=l ) (3-23) 

So the average velocity is 

v = (1 . 1 28)(7.785 x 1 0 3 )4f (3-24) 

The above is for the particles leaving the cathode, what we need, however, is 
the velocity distribution at the plane of the voltage minimum (the virtual cathode). 
It is calculated as follows. We know the energy in a particular bin is 2tijkT t , and 

for bin 14 for our temperature it is 1 .2808 eV. The factor of two occurs since we 
are considering those electrons that have escaped. The potential hill each particle 
must overcome is 0.655 eV (the voltage at the minimum is 0.655 volts). Thus, the 
fastest particles beyond the minimum have energy 1.2808 - .655 = .6258 eV. Now 
from experiments, we know the distribution is still Maxwell-Boltzmann, so this 
energy must correspond to those in bin 14 at the new, but lower, effective 
temperature. Bin 14 is characterized by rjj = 6.75. Thus our new temperature is 
defined by 

.6258 = (6.75 ]k Teff jkTeff 
Teff = 1 075" K 

Notice here the factor of 2 does not appear, as now we treat all emitted electrons 
as constituting a normal M-B gas. Thus our new temperature is 1075° K. We can 
always determine the effective temperature since the value of the voltage minimum 
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is known once the total current is chosen, and the saturation current is known for 
the specific cathode. The minimum voltage is given by 

V m = — (3-26) 

1 1,605 \ /J ) 

Where T is that for the cathode. Once the new temperature 7^ is known, the 
velocities for the bins become 

v. =3.893x10^7] T # ml sec (3-27) 

The new table for the velocities (at the voltage minimum) as compared to those 
in the Emission Modeling section (which were at the cathode), is: 


Bin 


v (m/sec) 

1 

.25 

3.191xl0 4 

2 

.75 

9.573x1 0 4 

3 

1.25 

1.596xl0 5 

4 

1.75 

2.234x1 0 5 

5 

2.25 

2.872x1 0 5 

6 

2.75 

3.5101xl0 5 

7 

3.25 

4.148xl0 5 

8 

3.75 

4.787x1 0 5 

9 

4.25 

5.425x1 0 5 

10 

4.75 

6.063x1 0 5 

11 

5.25 

6.7x1 0 5 

12 

5.75 

7.34x1 0 5 

13 

6.25 

7.98xl0 5 

_ 14 

6.75 

8.616xl0 5 
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To simplify the modeling, we may use a built-in feature of MAFIA; that 
being a parabolic distribution with chosen parameter. We will approximate the 
cosine distribution with the first two terms of its Taylor expansion, which is a 
parabolic representation. 


cos 6 = 1 + a 2 0 2 a 2 =- .4967 

;.p{d) = \-A9610 l 0<0< k / 2 


( 3 - 
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Beam Velocity Profile 

After a careful study of the experimental results on the beams created by a 
typical electron gun, we may make the following tentative statements. Firstly, 
the reported results are somewhat ambiguous, and secondly, they do not give the 
information in the form useful for our modeling effort. We desire the particle 
density and velocity at least over a specified plane normal to the axis of the beam. 
The reports give approximations that yield either velocity or charge density vs. 
radial value, but not both. Some studies have stated the current density falls off 
in a Gaussian manner from the beam center to the edge. Stating that the current 
density is Gaussian only tells us that the product of particle density and velocity 
are such. For our efforts, we need both quantities, not just their product. Some 
papers however, state that the beam is far from laminar (or Gaussian) and that rings 
of charge start at the edge and move through the beam to the center. It is not clear 
if the reported results are due to the particular electron gun used, or are typical of 
all guns. 

With the state of affairs as such, a reasonable approximation for the beam 
velocity profile chosen by C. Kory is as follows. At the emission plane initialize 
particles such that their angle with the normal 6 is given by 


tan# =— (1) 

U z 

Where u, and u- are the transverse and longitudinal velocities, respectively. The 
longitudinal velocity is constant, due solely to the anode potential. The transverse 
component is assumed to be Gaussian. Writing the above in more familiar 
notation 


y = tan 



) 


( 2 ) 
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In essence we are seeking the probability density function for 6 (or y), given the 
density for x (the transverse component u t ). 

When y > 0 

x = a tan y 


Let g(x) = y, then 


g'(x) 


a 


a 2 +x 2 


Then the density function for y (y = 6) is fy (y). 


a /{a + x~) a 


where fx(x) is the density of the random variable x (* = u,) 
If jc is normally distributed 




,S/2 c 2 


g42k 


Then 


fAy)= 


fl z +.v 2 1 

a g 42 k 


-a 2 tan 2 yllG 1 


When y < 0 * = a tan y 

So eq. (3) is the same as before. In eq. (3) 

2 2 i 2 

x =a tan y 


(3) 


Hence 


In most guns 


fy (y)=a- sec 2 y 


-a 2 tan 2 y /2<t“ 


<Tv27T 


(4) 
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u. 


u. 


io - 2 


\ — = 10 -2 = tany 

a 


;.y = e=.5iy 


The standard deviation is 


<7= — = 1.29x10 s V? 


m 


c = 


fr \ 

r 

( .24 } 

T 

95 ) 


•045 J 


= 28.44 


:.G = 6.88 xl0 5 m/ sec 

For typical values a ~ 2 x 10 s ml sec, a ~ 10 7 ml sec 

/ } .((9 = .573°)= 17.65 

Therefore eq. (4) is the probability density function for the emission angle for 
a given small patch, for a typical electron gun under the assumption of eq. (1). 
A sketch oifyiO) is below. 

a sec 2 0 


Me) 


-a 2 tan 2 0/2 ct 2 


<tV27T 


a-u_- 10 7 m/sec 



Figure 10 - The distribution f){0) for the emission angle <9 from the normal. 
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Conclusions 


It appears that MAFIA can be used to model electron guns with a Maxwell- 
Boltzmann energy distribution of the emitted electrons. To keep the run-time to a 
minimum, a scheme was devised that uses the same number of particles for any 
cathode current. The charge and mass of the particles are varied to conform to the 
given current. Similar to the older program EGUN, the calculation starts at a plane 
close to the actual cathode. In EGUN it is called the “starting surface” and is 
chosen somewhat arbitrarily (a few mesh cells from the cathode). EGUN sets up 
current tubes on the surface and varies them until the solution (forced to obey the 
Child-Langmuir 3/2 power law) is determined. We, however, emit charged 
particles (macroparticles) of varying charge and mass, which have the M-B energy 
distribution, and the proper distribution in emission angle from the cathode. Our 
emission plane is close to the virtual cathode whose position is known once the 
current is chosen. The actual virtual cathode varies with distance from the true 
cathode; it depends on the local cathode loading on a given patch. The loading 
varies by as much as 20% in some cases. The effective temperature of the 
macroparticles is calculated, given the known temperature of the cathode. The 
modeling of the beam should be closer to reality than any published scheme to 
date. The solution is not constrained to obey the 3/2 power law, so all regimes of 
emission are able to be examined. Specifically, if some region of the cathode is 
in saturation, this scheme can handle it easily. Also the non-laminar (noise) 
properties of the beam should be modeled completely. The use of particles will 
determine the true beam diameter as it exits the gun, as well as the spread of 
trajectories at the beam waist (near the point where the PPM stack applies the 
confining magnetic field). 
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The following bulletized list summarizes the accomplishments that are 
distributed throughout the memorandum. 

• Explained how the program EGUN operates, and some of its shortcomings. 
Fundamentally, the procedure bypasses the establishment of the space- 
charge layer, from which the anode current is developed. Instead, it models 
the beam as current tubes. The tubes cannot cross one another; so laminar 
flow is implied, even though experiments suggest otherwise. The noise 
properties of the beam are therefore lost completely. The calculation 
assumes the 3/2 law, which eliminates the possibility of portions of the 
cathode being in saturation. A consequence of the assumptions is its 
difficulty in correctly determining the beam diameter along the gun, as well 
as finding the correct value of the space-charge along the axis. An empirical 
factor of 5.5 is used to correct the value found by the calculation. It does, 
however, find the perveance to a few per cent. 

• Started with Lambert’s law as given in Gewartowski [6] 

j(6)=^cosd = qn(eXe)±q'£n j {e)v j (0) 

K 

and showed that n/0 ) is distributed as cosflfor a fixed value of velocity Vj. 
This result followed from the statement that the kinetic energy normal to the 
surface is twice that tangential to the surface. This was proven in Appendix 
I, and developed in Chapter 1, Emission Modeling. 

• With particle use, we will be able to model all parts of the I-V plane, as well 

as noise, and non-laminar flow. 1 - • ••• 

• Chapter 1 , Emission Modeling, indicates the basic way the problem should 
be handled. It ends by showing the numbers needed are too large for 
MAFIA. It demonstrated the difficulties encountered in trying to establish 
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the space-charge layer. It is actually this layer which serves to emit particles 
that form the beam. 

• It was recognized that the emission maybe viewed as follows. There is a 
M-B gas in the cathode, and there is another M-B gas, at a different 
temperature, which forms the space-charge layer. Then we can define a 
third temperature T e ff, which characterizes those electrons that have traversed 
the space charge layer and are going toward the anode. Those that are 
collected are also M-B which was determined from the experimental results 
of references [1], and [2]. Appendix I also showed the escaping electrons 

to have a temperature such that the average energy was 2kT vs the 3kT/2 of 
those that remain in the metal. Here T is the cathode temperature. However, 
in that calculation, the repelling force of the space-charge layer at the surface 
was not considered. 

• Chapter 2 detailed the solution for the 1 -dimensional diode that used the 
thermal velocity of the emitted charge as a boundary condition. The 
equations are non-linear, and tables have been established to facilitate the 
voltage, charge, electric field, and velocity field between the plates. We 
followed the lead of Beck [9] in this regard. 

• We stated the boundary conditions for the 3/2 law require the charge density 
to be infinite at the virtual cathode plane, which is not possible to be 
modeled using particles. 

• In principle, one should use particles in the following manner. With the 
cathode temperature given, the thermal limited current density is established 
for the particular cathode material. Emit particles with the M-B distribution, 
and which satisfy the current density constraint. The space-charge layer 
should form in the simulation, and steady state should occur. However, in 
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the middle range of the I-V characteristic, about Vi of those emitted will 
return to the cathode, only to be re-emitted. This would be a waste of 
computational time, just to establish this “bias condition” to allow the actual 
current to the anode to be established. Thus one must somehow sidestep the 
formation of the space-charge layer. 

• Chapter 3 outlined a method to sidestep the formation of the space-charge 
layer. The idea is to only emit those that reach the anode. The 
determination of their effective temperature T eff was presented. 

• Chapter 1 presented some bunches ( 1 7 angle sets) to be emitted if a 
distribution function in the program MAFIA was found to not be 
appropriate. A function f(0) was developed to determine the correct makeup 
of a given bunch. 

• We plan to model the cosine distribution for each energy bin group with a 
parabolic distribution function which is available in MAFIA. 

• Chapter 2 gave a solution for 1 mA of current and 1.172 volts. There, 
representative numbers were determined. In an actual gun, say the one in 
the Cassini Mission TWT, the current is 15 mA and the voltage is 5kV. 

Thus we must verify the model using the 1 . 1 72 case. 

• For the example we found the charge density varies by 4 orders of 
magnitude through the diode. The electric field changed from + to -, and 
varied by a factor of 22 from cathode to anode. The velocity field changed 
by a factor to 1000. Finally, the total number of charges between the plates 
in steady state is 166 million. Of these, 159 million are between the cathode 
and the plane of the voltage minimum. 

• Chapter 3 explained the method to keep the total number of particles in a 
given run constant. The variation in current is handled by adjusting the 
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charge and mass of the particles in the fixed 14 energy bins. This will 
dovetail into the program’s requirement of 4-5 macroparticles per mesh cell. 

• We used [1 1] to realize the average velocity near the voltage minimum is 
nearly constant for all current levels; it varies only slightly with cathode 
temperature. Looking back at the numerical solution obtained in Chapter 2 
we found complete compliance with this result! 

• With all of the compression done on the number of particles used, the runs 
may take 12-18 hours. We hope we have over-estimated the run times. 

• We have developed the probability distribution of the angle 6 from the beam 
axis for a reasonable approximation for defining an electron beam for 
other uses in parts of TWT simulations. If the transverse velocities are 
gaussian, while the normal component is constant, we find 6 is distributed as 


fr(Q)^fe(Q)=“ z SZC 2 d 


-uf /2<j 


<j42k 
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A ppendix I 

Section 5-15 of [8], page 141, concerns itself with the energies of escaping 
electrons from a hot cathode. This is an important result, as most formulas, results, 
concepts, etc. from electron gas theory in solids is for the case where the carriers 
remain in the metal. While we know the statistics of the ensemble for the electrons 
in the metal, what are they for the portion that escapes? It is more convenient to 
attack this problem as a function of velocity, thus 

r 1 2 
qE = —mv 
2 

Notice the notation that electric potential is given the symbol E (a peculiarity of 
Millman and Seely), while in Gewartowski it is given the label W. We use this 
notation since the development to follow is from Millman and Seely. Our goal is 
to show 


KE 


normal 


= 2 xKE 


tan gentiaJ 


( 1 - 1 ) 


The calculation proceeds by determining the average energy of those electrons in a 
metal that are capable of escaping. This is eq. (5-42) 

r r r I™( v =+v ! +v ! )w; 

E = - (1-2) 




After some calculations the following expression is found 



dv. 


Evaluating the integral in the numerator yields; let — = k t 

2 q 


then 


T k.vV h 'dx let u = v r 2 , du = 2vdv 

Jv B 1 x x xx 
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1 
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1 

t*T 
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B 

a 


2 A 


k\ •> 1 


A = 


2A 

m 


v“ H — 

VB X 


2qE T 

Their expression in the denominator is 

f v , e ~ Xv ' dv , 


d-3) 


and its evaluation proceeds as (a standard form) 


1 

~ 2X 6 


-E./Et 


E b = Min. of E 


The ratio is 


Now 


if we assume 


then 


“AVrfc-t -E B IE T 





m 

2 qE T 


V 


2 

xB 


E b 



2 

x B 


£~^b f£ T +E B IE t _j 


Then the ratio is 


(1-4) 


(1-5) 


d-6) 
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(1-7) 


. , u / 1 m , 1 m 2qE T _ _ 

*vi + Zx = T-v;» + - — — = E, + E r 

' A 2 q 2 q m 


as stated. The second term of eq. (5-42) is 


JJJ 


( 


1 m 2 Y 2 m 3 


v; \ —e E » ,ET Y^- Xv 'v x dv x dv y dv. 

[2^ y [ It r > 


111 


E w /£ r -Av;-AvJ-Av? tit 

— — e w r e “v av dv dv 

h 3 


X X y 


h^b,v;^ 


= ,, fc^L = *,-1- 

‘2A 




r ^ Xv> dv 

J-OO } 

/ 

2 


( 1 - 8 ) 


as stated. The final statements are: In leaving the metal, the average X-directed 
energy of each electron is reduced by an amount E B , so that the average X-directed 
energy of each electron after it has left the metal is ( E B + Ef) -E B = E T electron 
volts. Hence the total average energy of each electron that escapes from the metal 
is E t + 1/2 E t + \/2E t = 2E t electron volts. In other words the normal is E T 
whereas the Y-direeted (which is tangential here) is (l/2)E T . This is our first 
equation in this Appendix. 
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List of Symbols 


I d.c. gun (cathode) current 

K p perveance of electron gun 

J current density (amps/m 2 ) 

V voltage (or potential) along gun or 1 -dimensional diode 

E electric field in 1 -dimensional diode 

T \ cathode temperature 

T e ff effective temperature of the M-B ensemble that constitutes the anode 

current 

M-B Maxwell-Boltzmann energy distribution 

J 0 saturation current density at a given temperature for a specific cathode 

material 

0 angle from the normal to an emitting patch 

<J) azimuth angle 

v particle velocity 

< > average of quantity over the ensemble 

<v> average velocity of the particles over the ensemble 

W T thermal equivalent of energy 

m electron mass 

k Boltzmann’s constant 

q magnitude of electron charge 

n p (0) number of particles emitted vs. 0 (number/m 3 ) 

v(0) velocity of particles emitted vs. 0 (m/sec) 

U.S.A unit solid angle ( differential solid angle = 2 n sin0 d0) 
dP (u, 0) joint probability of emission velocity u and direction 0 
dN E number of particles with energy in the range E + dE 
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Njot total number of particles in ensemble (or a density :/m 3 ) 

E t W t energy (in Millman’s [8] notation). Note [8] uses E for both 

energy and voltage, a very confusing choice, even though they attempt 
to justify it in the beginning chapter 

p n dN n /drf Particle density with respect to the normalized energy r\ 

p E/E t normalized energy parameter 

E a constant = kT (Joules) or kT/q electron volts; for a temperature T 

Ipatch current from a patch of cathode 

Apatch area of patch 

nj particle density in energy binj 

ocj constant, j = 1,...,B 

B number of energy bins 

KE n0 nnai total kinetic energy normal to the surface 

KEtangentiai total kinetic energy parallel to the surface 

nj! particle density in binj with velocity at angle 0; 

K j a constant 

f(0) ycos 3 6 - cosd sin 2 6 

erf error function 


rij Ej/E t normalized energy parameter for binj 

r\ normalized parameter in the 1 -dimensional diode; not to be confused with 

r| for normalized energy 
r\ c T) parameter at cathode 

£ normalized distance parameter defined by Langmuir [10] 

^c <a normalized parameter at cathode or anode 

d separation of planes in 1 -dimensional diode 
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V m 

voltage at minimum (virtual cathode plane) 

V a 

voltage at anode 

K 

constant in analytical expression for E 

PW 

charge density (Coulombs/m 3 ) 

£ 0 

permittivity of free space 

A 

cross-sectional area 

x a 

position of anode 

v(x) 

velocity field in the 1 -dimensional diode 

TWT 

traveling wave tube 

tf 

time of flight of a particle in the 3/2-law regime 

ip 

patch current 

Np 

number of patches 

AQ P 

total charge emitted from a patch during At 

AQj 

net charge in binj 

n p 

number of particles/patch emitted during At 

q» 

charge of particles in binj 

Mj 

mass of particles in binj 

m 

mass (generally an electron) 

n j 

number of particles (per m 3 ) in binj 

u t 

transverse velocity in electron beam 

Uz 

longitudinal velocity in electron beam 

EGUN 

popular electron gun program (used for last 25 years) 

PPM 

periodic permanent magnet (stack) 
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